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Abstract 

Interacting particle systems are continuous time Markov processes 
which are used to construct stochastic spatial models. Monotonicity is 
a useful property which simplifies certain calculations, one of which is 
the ability to use computational algorithms to sample exactly from the 
stationary distribution for certain processes. A monotone interacting par- 
ticle system is called attractive. Monotonicity is well understood for spin 
systems which only include two particle types, such as the contact pro- 
cess, however, when constructing applied models, it is often desirable to 
include more. In this paper, an interaction map is used to describe the 
interactions that occur in a model and to understand monotonicity for a 
certain class of multitype contact processes. 

1 Introduction 

Interacting particle systems are a class of Markov processes used to model the 
evolution of particles types S on a collection of sites A C 7L d . Consequently, a 
state for these processes is a mapping r\ : A — > S that assigns a particle type 
to each site. A process realization r] t = {rj t : t e [0, oo)} is a right continuous 
function oft. Elements of the state space, S = S A , are interchangeably referred 
to as states or configurations. If 7](x) = a, site x is said to be infected or 
occupied by particle type a, or more simply, particle a sits at site x. Typically 
£ 5 and r](x) — indicates an unoccupied site. 

A general interacting particle system is quite difficult to analyze and some 
additional properties must be considered in order to make meaningful general 
statements. In this paper, our motivating examples are multi-species biological 
models. Thus, we will be primarily focused on circumstances in which the 
number of particles types #\S\ is at least three - vacant and two species of 
organism. 
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The contact process [5] is a basic interacting particle system used to build 
many biological models, and the nature of the interactions which occur has led us 
to introduce the 'interaction map' which characterizes the allowable interactions 
in a model. 

The main result is a characterization of those interaction maps that lead to 
monotone processes. In this case, the particle types S must be ordered thus 
endowing the state space with a partial ordering '<' (sitewise comparison). 
Monotonicity allows us to reduce the analysis of many properties of an interact- 
ing particle system to a study of the process considering the evolution starting 
from the small number on initial configurations that are extremal under the 
given partial order. For example, the coupling from the past (CFTP) algorithm 
[TT] allows us to obtain an exact sample from the stationary distribution of an 
irreducible and aperiodic Markov chain with a finite state space. For monotone 
processes, the algorithm need only be applied to the extremal states and n (the 
states where all sites are identically occupied by or n particles respectively) 
such that < r] < n for all r\ e S. 

The interaction map formulation also allows for a fast assessment of whether 
or not reordering the particles may make the process attractive. We show that 
switching the order of the basic two type contact process results in attractive- 
ness. 

We will briefly review the contact process and monotonicity. This will set 
the stage for the development of the interaction map. 

2 Contact Processes 

The contact process is a spin system (it only has two particle types), and it is 
the basic interacting particle system of interest here. Its state space is {0, 1} Z 
with particle type representing an empty site and type 1 an occupied site. 
The transition rates are: 



where ri\ is the number of occupied neighbors. A birth occurs at a rate propor- 
tional to the number of infected neighbors. An occupied site becomes vacant 
at rate one. The transition rate for a spin system is referred to as the flip rate, 
and denoted c(x,ij). Thus for the contact process, 



where Ii(-) is the indicator function for particle type 1. 

Numerous variations have been created around the rule that infection occurs 
at a rate proportional to the number of occupied neighboring sites. A contact 
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process which has n + 1 particle types will be referred to as an n-type contact 
process (usually type denotes an unoccupied site). The basic two- type contact 
process is a competition model studied by Neuhauser [10] . This process has the 
following transition rules: 

— > 1 Ai^i — > 2 X 2 n 2 
l->0 Si 2^0 S 2 

The two-stage contact process studied by Krone [7j is another two-type contact 
process. It is a single species model which includes two life stages. The transi- 
tions are such that 2's give birth to l's and the l's mature into 2's at a constant 
rate: 

0-^1 \n 2 1 -» 2 7 
U0 1 + 5 2^0 1 

The grass-bushes-trees successional model proposed by Durrett and Swindle [2] 
is a two-type contact process as well. It is the basic two-type contact process 
with the modification that 2's are allow to give birth onto sites occupied by l's 
in addition to empty sites: 

— s- 1 Aini — > 2 A 2 n 2 1 — > 2 A 2 n 2 

1 ->• Si 2^0 S 2 

These models are all multitype contact processes, however they are no longer 
spin systems so what is generally required for attractiveness has not been pre- 
viously known. In [7] (Theorem 3.1), it is shown that the two-stage contact 
process is monotone with respect to its parameters. In [T] (Proposition 1.1), a 
more complicated type of monotonicity property is described for a three-type 
contact process. 



3 Monotonicity 

A Markov process with a partially ordered state space, S, and semigroup, S(t), 



is called monotone if either of the equivalent conditions, (2a I or (2b), is shown 
to be satisfied. 

/ G J( implies S(t)f € M for all t > (2a) 
Hi < fj, 2 implies fiiS(t) < fi 2 S{t) for all t > (2b) 

The fii are probability measures on the state space, and ^# is the set of con- 
tinuous monotone functions, / : S — > M, such that two states satisfying r\ < £ 
implies f(rj) < /(£). For the proof of equivalency, see [8] (Chapter 2, Theorem 
2.2). Furthermore, fii < \i 2 is equivalent to there existing a measure, 1/onSxS, 
that satisfies 



(a) z/{(?7,0 :ijeA}=ft 1 (A), and 

(b) u{(tj, £) : ^ G A} = [12(A), where A is any Borel set in S, and 
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(c) f{fa,0 :»?<£} = !• 

For the proof of this, see [5] (Chapter 2, Theorem 2.4). For interacting particle 
systems, the typical route is to define what it means for the process to be 
attractive and then to prove that this is equivalent to being monotone. 

For a spin system, given 77 < £, attractiveness is equivalent the the flip rate 
satisfying 

c(x, T)) < c(x, £) when r)(x) = £(x) = 0, (3a) 
c{x,rj) > c(a;,£) when tj(x) — £(x) = 1. (3b) 



See [5] (Chapter 3, Theorem 2.2) for the short proof that (3a) and (3b I are 
together equivalent to monotonicity. This paper generalizes these conditions to 
a class of multitype contact processes. 



3.1 The Interaction Map 

The approach here is to define what is called the interaction map. We start 
with a finite set of totally ordered particle types, S = {0, 1, 2, . . . , n} and define 
a mapping that describes all interactions between them. 

Definition 1. Given a finite set of totally ordered particle types, S, the map, 
J" : S x S — > S is called an interaction map if its domain is all of S x S, and its 
range is a subset of S. The particle type that replaces type a upon interaction 
with type b is given by J* (a, b), and this interaction is denoted by a \ J? (a, b). 

The domain of J! is partitioned into sets of up, null, and down interactions, 
respectively: U = {(a,b) eSx S\J(a,b) > a}, N = {{a,b) eSx S\j?(a,b) = 
a}, and V = {(a,b) £ S x S\ c p(a, b) < a}. A process whose interactions are 
completely defined by a single interaction map is referred to as a interaction map 
system (IMS). Many processes may also be described using multiple interaction 
maps; this is discussed briefly in the last section but is not pursued in detail 
here. 

For the contact process, ^(0, 1) = 1 denotes a birth onto an unoccupied site 

1 i 

(0 — > 1), and ^{\,i) = for any i denotes a death (1 — > 0). The interaction 
map for the voter model is the same as for the contact process except that 
^(1, 1) = 1 because a voter only changes opinions by contact with its opposite 
(Figure [lj . The interaction map for the basic two-type contact process is built 
from that for the contact process. It additionally has ^ (0, 2) = 2 and ^ (2, i) — 
(for any i) for births and deaths of species two respectively (Figure [2]). 

The interaction map, J! , is called non- decreasing on A C S x S if for any 
(di,&i) and (02,62) hi A which satisfy (a\,b\) < (02,^2) (meaning a\ < and 
h < b 2 ), it follows that J r (a 1 ,b 1 ) < ^(a 2 ,& 2 )- 

Definition 2. The interaction map will be called attractive if it satisfies the 
following conditions: 

(a) Jf is non-decreasing on U. 
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(c) (v) 
Figure 1: Interaction maps for the (c) Contact Process and (v) Voter model. 








1 


2 














1 


1 








2 


2 









Figure 2: Interaction map formulation for the multitype contact process 



(b) ^ is non-decreasing on V. 

(c) If (ai,&i) < (a 2 ,6 2 ), (a>i,h) eU, and (o 2 ,6 2 ) £ VUAf, then J {a Xl b x ) < 
a 2 . 

(d) If (ai,6i) < (02,62), (ai,fei) £ U U Af, and (02,62) € £>, then a x < 
(o 2 ,6 2 ). 



The last two conditions of Definition [2] are equivalent to: 
(c 1 ) If (01,61) £ W, then /{axM)< min{o : (a, 6) eDUiV}. 
(d') If (a 2 ,6 2 ) € X>, then max{a : (a, 6) e U U TV} < J (a^bi). 

a 

Definition [2] is designed so that given two ordered pairs of particle types, 
(01,61) < (02,62), we have ^"(01,61) < ^(02,62) when both pairs have an 
up or both a down interaction (or one of them is null). When one interaction 
is up and the other is down, the ordering need not be preserved, but cannot 
be broken arbitrarily; it must obey (c) and (d) of Definition [2] (or equivalently 
(c') and (d')). Staring from Definition [2] a little extra work shows that an 
attractive interaction map ^(0,6) is nondecreasing in 6 and is nondecreasing 
in a except for particle pairs satisfying (a, 6) G U and (a + 1, 6) G T>, in which 
case ^ (a, b) = a + I and ^(a + 1, 6) = a. 

The interaction map of the basic two-type contact process (Figure [2]) is not 
attractive because /(0, 2) = 2 > min{a : (a, 2) £ AfWD} = 1. There are several 

a 

ways to modify this interaction map into one which is attractive. Keeping the 
birth onto empty sites unchanged, necessitates that (a) ^(1,2) = I or (b) 
^(1,2) = 2 and ^(2,2) £ {1,2} (Figure [3). Thus the death rates are no 
longer both constant, or the two species competition character of the model 
may be broken. 
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(a) (b) (c) 

Figure 3: Attractive modifications to the basic two-type contact process inter- 
action map. The interaction between (1,2) and (2,2) must be modified so that 
Definition [2] is satisfied. The death rates are no longer constant. 



Figure[4]shows four examples of attractive interaction maps with the different 
set of interactions shaded. Dark gray denotes down interactions, light gray 
denotes up interactions, and no shading denotes null interactions. We can see 
that the interaction maps are nondecreasing separately in the light gray and 
dark gray regions and that when a light gray cell is on the left of a dark gray 
cell, the decrease is exactly one: (a,b) = a + 1 and J* (a + 1, b) = a. 
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Figure 4: Examples of attractive interaction maps are shown above. Those cells 
shaded dark gray represent down interactions, light gray up interactions, and 
unshaded cells are null interactions. For example, in the upper left interaction 
map, ^(2, 3) = 3. This shows that particle type 3 influences type 2 to become 
type 3 and is thus an up interaction. 



4 Transition Rates 

The transition rates are denoted by r u (r],x,y) and rd(ri,x,y), which give the 
rates for up and down interactions respectively. By convention, null interactions 
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are assigned a transition rate of zero since they do not change the configuration. 
For each pair of particle types (a, b) e S x S, \ a b denotes the rate at which the 

interaction a \ J* {a,b) occurs. For a configuration 77 the rates at which a 
particle at site y influences a particle at site x to change are thus given by 

r u (i],x,y) = *ab<l>(x,y)I{(a,b)}(v(x),v{y)) 
aes,beu a 

rd(v, x ,y) = X! ^ab<t>(x,y)I{(a,b)}(v( x ),V(y))- 

aeS,be"D a 

The function cf>(x,y) defines the neighborhood, and is assumed to be non- 
negative and bounded over a finite neighborhood. The sets U a = {b € S : 
^(a,b) > a} and T> a = {b E S : J" (a,b) < a} for a given particle type a 
denote the sets of particle types with up or down interactions respectively. 

The total transition rate will be denoted r(r), x, y) = r u (r), x, y) + r<f (77, x, y). 
Note that r u and are never simultaneously nonzero for any given configuration 
since we are restricting ourselves to processes which are described by a single 
interaction map. This is not true in general for multitype contact processes. 

The generator for an IMS is defined on / 6 C(S) by 

om = E r ^ x > y)(f(v xy ) - /(»?))■ (5) 

The state rj xy represents the state 77 with the particle at site x changed according 
to the influence of the particle at site y: 

{fi(z) if z^i 

/{r){x),r]{y)) if z = x 

4.1 The Coupled Rates 

Suppose we have two interaction map systems, r] t and £t on the state space 
S, with transition rates r\(r),x,y) and r 2 {^,x,y) respectively, with ri{-,x.y) — 
r iu{',x,y) + rid(-,x,y). The coupling described below is essentially the same as 
the Vasershtein coupling [T^] given for spin systems, also known as the basic 
coupling [8]. The coupled process, (77t,£t), is a Feller process whose state space 
is 5 x 5 and evolves according to the following rates: 

!(ri xy ,£ xv ) at rate r(r),£,x,y) 
(r),£ xy ) at rate r 2 (£,x,y) - r{r),£,x,y) (6) 
(v xy ,0 at rate r x (77, x, y) - r (77, £, x, y) 
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where f (77, £, i, y) = min(r lu (r], x, y), r 2u (f , x, y)) + min(ri d (?7, x, y), r 2d (£, x, y)). 
The generator for the coupled process is defined for / G C(S x H) by 

Gf(v,0= E iri(v,x,v)-Hv,Z,x,v))V(rr,Q-f(v,£)) (7a) 
+ E ( r ^,x, y )-r( V ,^x,y))(f( V ,e v )-f(v,0) (7b) 

+ E Kv,Z,x,y)(f(ri xy ,t xv ) -f(v, 0) (7c) 

This coupling is shown to preserve the partial ordering under certain conditions. 



5 Monotonicity and Attractiveness 

Now we are ready to prove that the coupling preserves the partial ordering 
on the state space almost surely when the interaction map is attractive and 
transition rates satisfy certain conditions. This next proof is almost identical to 
that of Theorem 1.5 in Chapter 3 of [5], but since the processes here are slightly 
more complicated a bit more work is necessary. 

Theorem 3. (Extension of Theorem 1.5 in Chapter 3 of [5]) Suppose that the 
processes 7?t and £t on the state space S share the same attractive interaction 
map. Furthermore suppose that whenever 77 < £, the transition rates satisfy 

riu(v,x,y) < r 2u {£,,x,y) when rj xy (x) > £(x), (8a) 
rid(v,x,y)>r 2d (^x,y) when ^ y (x) < T)(x). (8b) 

Then for all 77 < f and t > 0, 

P {r > 4) fa = (9) 

Proof. Suppose we have states 77 < £, then all that is necessary is to show that 
the coupled process preserves the partial ordering almost surely. Note that the 
coupling only allows simultaneous transitions if they are both up or both down. 



Case 1 If the lower configuration may jump above the upper, £(x) < rj xv (x), 



then r lu (rj,x,y) < r 2u (£,x,y) by (8a). Since 



r(v> £> x, y) = min(ri„(?7, x, y), r 2 „(£, x, y)), 

we get 7*1(77, x, y) — r(r), £, x, y) = 0. So the problem transition, (77, £) — > (r] xy , £), 
occurs at rate zero. 



Case 2 If the upper configuration may jump below the lower, £ xy (x) < r](x), 
then ru(ri,x,y) > r 2d (£,x,y) by fib] ). Then r 2 (ij,x,y) -r(rj,£,x,y) = since 
r 2 is the minimum here leading to (77, £) — > (77, £ xv ) at rate zero. 
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Case 3 If £(x) < £ XJ/ (x) and rj(x) < r) xy (x), then rf v (x) < £ xy (x) by the first 
attractive interaction map property. Assuming that £(x) < rf v {x) would put 
us back in case 1 and we are done. If r] xv (x) < £(x), then there are no problem 
transitions that break the partial ordering. 

Case 4 If ^(x) < £(x) and rf y (x) < T)(x), then rf v {x) < £ xy (x) by the 
second attractive interaction map property. If ^ xy (x) < rj{x), then we are back 
in case 2, otherwise all transitions preserve the partial order. 

Case 5 If r](x) < rf y (x) and £, xy (x) < £(x), then ?f y (x) < £(x) by the third 
attractive interaction map property. So r(T),£,x,y) — since the minimum 
up transition rate is as is the minimum down transition rate. This gives 
(77, £) — > (rj xy ,^ xy ) at rate zero, and the order is preserved. 

Case 6 If r/(x) < rf y (x) and ^(x) < £(x), then r)(x) < £ xy (x) by the fourth 
attractive interaction map property. So f(j], £, x, y) = since the minimum up 
transition rate is as is the minimum down transition rate. This gives and 
(77, £) — > {rf y ^ xy ) at rate zero once again preserving the partial ordering. 

□ 

Definition 4. An IMS will be called attractive if it has an attractive particle 
interaction map and given rj < £, the transition rates satisfy: 

r u (v, x, y) < r u (€, x, y) when J{r\{y), r/(x)) > £(x), 
rd(v, x , y) > rd(€, x, y) when J{i{y), £(x)) < ry(x). 

Consider two ordered configurations, 77 < ^. In the event that an interaction 
may cause 77 to jump above £ at the site x, the latter configuration must have an 
equal or larger up transition rate and a jump transition that goes at least as far: 
£(x) < ri xy {x) < £ xy (x). When the upper configuration, £, could possibly jump 
below the lower configuration, 77, a similar statement applies. A coupling which 
preserves the partial ordering of the underlying processes is called a monotone 
coupling. 

Theorem 5. (Extension of Theorem 2.2 in Chapter 3 of [H]): An IMS is mono- 
tone if and only if it is attractive. 

Proof. If the process is attractive, then we have a coupling which preserves the 
partial order that proves that the IMS is monotone. Assuming that the process 
is monotone we must prove that it is attractive. We fix x and choose 77 and £ 
such that r](y) = for all y ^ x. Choose b such that -q(x) < £(x) < b. This 
will be used to prove conditions on up transitions. If no such b exists, there is 
no problem, as T)(x) cannot jump above £(x). Similarly if there is a b such that 
b < rj(x) < £(x), then this is used to show something about down transitions. 
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We start with the monotone function fbx(v) — I[b..n](v( x ))> the indicator on 
all particle values bigger than or equal to b. Since fb x is monotone, S(t)fb x (v) < 
S(t)fbx(0- Noting that b was chosen so that fbx(v) = fbx(0 shows that: 

s(t)fb X (v)-f bx (v) < s(t)f bx (Q-f bx (0 
t ~ t 

Taking the limit t — > gives: 

Gf bx (v)<Gf bx (0- 
Plugging in the form of the generator [5] results in: 

r ( r t^,y)(f bx (v xv )- f bx {v)) 

y£N(x) 

< ]T r(tx,y)(f bx (e v )-fbx(0) 

y€N( x ) 



(11) 



This gives the total rate that rj(x) transitions up into the set [b..n] = {a£ S\b < 
a < n} is less than or equal to the total rate that £(:r) goes up into the set [b..n], 
and the total rate that rj(x) leaves the set [b..n] is greater than or equal to the 
total rate that £(x) does the same. 



Then, due to the choice of transition rates, r(r), x, y) = \ v ( x ) tV ( y -)(j)(x, y), (11 1 
becomes 

(12) 



< 



Because the same particle sits at all y, J2 y 4>(x, y) cancels out from both sides. 

\(x)Mv)(fbx(v xy ) ~ fbxiv)) , 13 s 



The following inferences come from looking at the possibilities for ( 13 ) 

(i) If rj(x) < £(x) < b, then 

(a) rf«{x) >b^£ xy (x) > b 

(b) £ xy (x) <b^rf y {x) <b 

(ii) If b < r)(x) < then 

(a) C xy (x) <b^rf y {x) < b 

(b) rj xy (x) >b^£ xy (x) > b 



Letting b = £(x) + 1 in inference (i,a), we see that when r) xy (x) > ( [13] 
becomes \j(x),ri(y) < ^Z(x)£(y) which is our first rate restriction. For the seconc 
rate restriction, let b = rj{x) in inference (ii,a), we see that when rj(x) > £ xy (x) 
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(13) becomes X v (x).ri(y) > ^Z(x),£(v) finishing the rate restrictions for being an 
attractive process. 

If both states have possible up transitions, r](x) < rf v {x) and £(x) < £ xy (x), 
then letting b = rf v (x) shows that under (i,a), rf v {x) < £ xy {x). Similarly if 
both states have possible down transitions, rj xy (x) < rj(x) and £ xy (x) < £(x), 
then b = £ xy (x) + 1 under (ii,a) shows that rf v (x) < £ xy (x). These give the 
first two requirements for an attractive interaction map. 

Assuming rj{x) < t] xy (x), £ xy (x) < £(x), and b = £(x) + 1 along with (i,b) 
shows that rf y {x) < b. If rj(x) < rf y {x) and £ xy {x) <£(x), then letting b = r](x) 
shows that b < £ xy (x) by (ii,b). Now the last two attractive interaction map 
requirements are met. Since this did not depend on the particular choices of 
particle types that sat at x and y, we are done. □ 



The proof of Theorem [5] is somewhat more involved than the case for spin 
systems since particle values are allowed to increase by more than one in an 
interaction giving rise to the possibility of 'jump overs'. This result applies to 
any interaction map system. The inclusion of spatial and temporal variations 
in <j)(x, y) or in the rate parameters does not affect this result so long as the 
conditions for being attractive hold over the entire lattice at all times. 



5.1 Reordering the Particles 

If one develops a model which either does not have an attractive interaction 
map, or the rate restrictions which allow this model to be attractive are not 
desirable, a re-ordering of the particle values may give an attractive model or 
more desirable rate restrictions. 

The basic two-type contact process can be made attractive with a particle 
re-ordering. The issue that the interaction map is not attractive is resolved by 
making the permutation {0, 1, 2} — > {1, 0, 2}, in the sense that now 1 < < 2, 
the interaction map is then attractive. To avoid confusion on such an awkward 
ordering of integers, we label particle as one species (species 0), particle 1 as 
empty, and particle 2 remains labeled species 2. The transition rates are still 
as before: empty sites become species i at rate /3j times the number of species i 
nearby, and species i dies at constant rate Si. This information is summarized 
in Figure [5] 
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Figure 5: Reordering of the multitype contact process for attractiveness. 
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This model now has an attractive interaction map and no rate restrictions. 
Normally represents an empty site, but this shows that thinking more carefully 
about particle labels is beneficial. This gives us ordered extremal stationary 
distributions: < S 1 < v 2 where 8 l is point mass on the state rj(x) = i for all 
x and Vi = lim 5 l S(t) is the invariant measure for species i in isolation. 

t— f oo 

6 The Graphical Representation 

Now we define a graphical coupling which couples all states simultaneously. This 
graphical representation was first introduced by Harris [U [6] and subsequently 
modified by others [3} [HI IS] - This method is useful for exact sampling algorithms 
such as CFTP. THe graphical coupling here is built from that in [3J. 

For each site, x, and every y such that cj>(x,y) > 0, let U xy and D xy be two 
independent, identically distributed Poisson point processes on (0, oo) x (0, oo) 
with intensity equal to two-dimensional Lebesgue measure. Assume that our 
rates have an upper bound, c. For each x, define 1 X = {T x i < T x 2 < ■■■} by 
T Xi o — and 

T x<n = inf{f > T X) „_! : (v, t) £ \J U xy U D xy for some v < c} 

y 

which will represent the times at which transitions could possibly occur at the 
site x. This is a projection of a union of independent Poisson point processes, 
and the intensity measure of points in % x is 2c=ffN(x, p) where N(x,p) = {y : 
\x — y\ < p and cf)(x, y) > 0}, so 1 X is also a Poisson point process. 

The graphical representation is created on a grid A x [0, oo). At each point, 
t, in if 3u,y such that (it, t) £ U xy draw an arrow from y pointing to x 
with an open circle at x. If 3u,y such that (u, t) € D xy draw an arrow from 
y pointing to x with a closed circle at x. Write the u values next to the tip of 
each arrow. Figure [6] give a realization of the graphical representation for a one 
dimensional index set and rates bound above by c = 3. 

These point processes are used to evolve the interaction map system. Given 
an initial state, 77, the following theorem describes how the graphical coupling 
evolves the process over time. 

Theorem 6. The path r)t is constructed according to the following rules is the 
interaction map system with interaction map ^ , and generator given by ([5]). 

(a) Up transition rule: The particle at site, x, is replaced by the particle 
given by J ' (rj t -{x) 1 r]t-{y)) at time t £ 1 X if there exists a u such that 
(u,t) £ U xy , with c /(r)t-(x),rit-(y)) > r)t-(x), and u < r u (r] t _, x, y). 

(b) Down transition rule: The particle at site, x, is replaced by the particle 
given by J! [rj t _{x),ri t _{y)) at time t £ £ x , if there exists a u such that 
(u,t) £ D xy , with c /(r]t-(x),rit-(y)) < r) t -(x), and u < r d (rj t ^,x 1 y). 
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Figure 6: A realization of the graphical representation of the point process 
coupling. All points with u > 3 have been left out. Closed and open circles 
represent possible down and up transition points respectively. 



The construction of these Poisson point processes and Theorem [6] is based 
upon the construction in Chapter 32 of [3]. This is the basis for an accurate 
method of simulating these types of processes. 

Theorem 7. The graphical construction in Theorem [6] maintains the partial 
order of an attractive IMS. 

Theorem [7] is proven by the fact that all transitions preserve the order of 
configurations. 

Theorem 8. Consider two sets of parameters satisfying }P) < A^? for all up 

transitions and > A for all down transitions. If two states satisfying £W < 
£( 2 ) are the initial states for the processes with the corresponding parameter sets 
above, then £^ < £^ for alH > for the above graphical construction. 

Theorem [8] is proven by applying Theorem [3] to the graphical coupling. This 
shows that an attractive interaction map process is monotone in each of its rate 
parameters. This is similar to the monotonicity theorems in [T] and [7J. 
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7 Multiple Interaction Maps 



One may desire to include more complicated interactions such as non-unique 
interaction, i.e. when a particular particle type can have multiple influences 
over another, or if a constant transition rate is included along with several other 
particle interactions. Take for example the grass-bushes-trees model studied in 

In order to make this model attractive, we reorder the particles as for the 
basic two-type contact process and introduce a separate interaction map for the 
births of trees. The basic interaction map, ^i, and transition rates are given by 
Figure [7j This interaction map is the basic two- type contact process map with 
births for species two removed and is attractive with no rate restrictions. Now 
we just need to account for this extra birth event. This amounts to including 
the extra interaction map, given by Figure [HJ and this map is attractive 
as well. This is not a unique formulation, deaths could be distributed among 
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Figure 7: Grass-bushes-trees initial rate parameters and interaction map, 
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Figure 8: Extra rate parameters and interaction map, for the birth of trees. 

both interaction maps and still maintain attractiveness. When simulating this 
process we need three Poisson Point Processes: U = up transitions for which 
includes deaths for species zero, D = down transitions for ^\ which includes 
deaths of species two and births of species zero, and B = birth events for species 
two (up transitions for ^2). This shows that the grass-bushes-trees model is 
attractive with no rate restrictions. 

The idea is that any number of interaction maps can be used. If each map is 
attractive, then a collection of rate restrictions allows the model to be monotone. 
Each map is assigned distinct up and down point processes associated with it 
in the graphical coupling. The equivalency of attractiveness and monotonicity 
for a process not described by a single interaction map is not discussed here in 
detail, but a sufficiency condition is given. 
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Theorem 9. Monotonicity Sufficiency Condition: If an interacting par- 
ticle system is formulated with multiple attractive interaction maps, J"i i — 
1, ...,n, and the corresponding attractive transition rates, then {{7j, -Di}i=i,...,n 
is a monotone coupling for the process where {/, and -D^ are the point processes 
for the up and down transitions for interaction map respectively. 

Each of the point processes \U\,D\,XJi,'Di, —,U n , D n } preserves the partial 
ordering of the state space, thus we have a monotone coupling. This allows non- 
unique interactions between two particle types. If multiple interaction maps are 
not used, monotonicity is given by inequalities involving sums of rate param- 
eters rather than a comparison of individual parameters. This may allow the 
relaxation of the requirement that each individual map and its parameters be 
attractive, but this is not pursued further here. The grass-bushes-trees model 
in Figures [7] and [8] is attractive with no rate parameter restrictions according to 
this theorem. The two-stage contact process can also be seen to be monotone 
with no rate restrictions. 

With the interaction map formulation presented here, monotone properties 
of multitype contact processes can be assessed quickly. While the main result 
only applies to processes with a single interaction map, it is still useful for de- 
termining whether or not a process with multiple interaction maps is monotone. 
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